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Abstract: We split the Yang-Mills Lagrangian into a free and an interaction part in 
such a way, that the free part is non-local and contains an arbitrary form factor. Manifest 
gauge invariance is guaranteed by connecting the field-strength tensors at different space- 
time points by a string. As a result the gluon propagator, which, due to the presence 
of the string, now contains many different contributions, comes out strictly transversal 
in one-loop order. We also calculate the ghost self-energy and the ghost-gluon vertex 
in one-loop order. Subsequently we discuss how one can determine the "optimal" form 
of the form factor. To apply well known principles like "fastest apparent convergence" 
or "principle of minimal sensitivity" one has first to solve some problems connected 
with divergences and renormalization. Here we concentrate on the calculation of the 
anomalous dimensions and the /^-function. This is technically simpler because only the 
divergent contributions of the integrals have to be determined. The /3-function becomes 
independent of the gauge parameter as it should. A puzzle with respect to the principle 
of minimal sensitivity shows up. Really interesting non-perturbative results are expected 
when applying the above principles directly to the propagator. For the /^-function a two 
loop calculation would be required to obtain non-trivial results. 
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1 Introduction 



The "optimized ^-expansion" , also called "linear ^-expansion" or, more appropriately, 
"variation perturbation theory", is a powerful method which combines the merits of 
perturbation theory with those of variational approaches. The underlying idea is simple. 
Generically, the Lagrangian is split into a free and an interacting part in such a way 
that an arbitrary parameter A (or more) is artificially introduced. The interacting part 
is multiplied by a factor 5 which serves as formal expansion parameter and is put equal 
to one at the end. The exact solution should be independent of the parameter A while 
any approximate solution will, of course, depend on it. One way to fix the value of 
A is the "principle of minimal sensitivity" (PMS): It demands that the approximate 
solution should depend as little as possible on the parameter. This means that A should 
be chosen such that the quantity to be calculated has an extremum. In this way the 
result becomes non-perturbative because A becomes a non-linear function of the coupling 
constant. In every order of perturbation theory the optimal value of the parameter A 
has to be calculated again. 

The method is now well established and it is therefore impossible to give all references. 
For the older literature we refer to Stevenson [jlj and references therein, some more recent 
references can be found, e.g. in [0. 

The field-theoretical applications have been rather modest up to now as far as the 
splitting of the Lagrangian into a free and an interacting part is concerned. Essen- 
tially only the mass parameter was used as variational parameter in the free part of 
the Lagrangian. The method is then also called "Gaussian effective potential". Some 
references are given in ||. In 0, on the other hand, the whole Lagrangian was scaled 
with a constant. 

It has been suggested in the literature || that one should start the procedure with the 
most general free Lagrangian. This would lead to a non-local quadratic action containing 
an arbitrary form factor. It is also believed that the behavior in the ultraviolet region 
should be identical to that of usual perturbation theory if the form factor approaches 
one for large momenta. We will see that the situation is in fact more complicated in our 
case. 

In the present paper we apply this idea to QCD, for simplicity without quarks for the 
present. The first central point, treated in sect. 2, is the construction of a Lagrangian 
which is gauge invariant for every 5. This can be done by connecting the field tensors, 
now taken apart to different space time points in the non local Lagrangian, by a path 
ordered string. Gauge fixing can be done in analogy to the familiar case. The action 
becomes an infinite series in 5 which coincides with the usual one for 5 = 1. For explicit 
calculations one expands the action up to the desired order. 

In sect. 3 we give the general formula for the gluonic vacuum polarization. The 
transversality can be explicitely checked, thus confirming the manifest gauge invariance 
of the formalism. In sect. 4 we give the corresponding expressions for the ghost self- 
energy and the ghost-gluon vertex. 

Up to this point the approach is quite general and the form factor, introduced when 
splitting the Lagrangian, essentially arbitrary. In sect 5. we discuss some conceptional 
questions which arise when one tries to fix the form factor in an optimal way. Due to 
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divergences this is a non-trivial task. In sect 6. we therefore specialize to a scale-invariant 
ansatz and concentrate on a calculation of the anomalous dimensions and the /3-function. 
This is much simpler because only the divergent contributions have to be extracted. 

Due to the complexity of the problem we have not yet achieved a break-through in 
the present paper. We believe, however, that we have demonstrated the feasibility and 
the potential of the method. With some optimism one may hope that it can become a 
new alternative approach to quantum field theories. In the conclusion we discuss, how 
non-perturbative results could be obtained in future applications. 

2 The action 

Our starting point is the classical Yang-Mills action 

S {d) = ~\J F^(x)F' M/a {x)dx. (2.1) 

The only modification at this point is that we replace the bare coupling constant go in F£ v 
by 5 go. The same replacement is made in the definition of gauge transformations. This 
is necessary in order to obtain gauge invariance for any 8, which is of central importance. 
Thus we put 

F; u = d,Al - d u Al + bg«r hc A\Al = G% + 5g Q C^. (2.2) 

For later convenience we have introduced the notations G a ^ u and C^ u for the abelian part 
and the commutator term of respectively. 

In order to start the 5-expansion with a general non-local free action we take the two 
factors F® u (x) in (2.1) apart to different space time points x,y and introduce a Lorentz 
invariant form factor K[x — y). To do this in a gauge invariant way one has to connect 
the points x, y by a path ordered octet exponential U(x,y): 

U(x,y) = Pexp[% [' V" A«(z)dz% (2.3) 

with (T a ) bc = —if abc the representation matrices in the octet representation. Tech- 
nically this is much simpler than using two factors U (x, y) and U (y, x) in the triplet 
representation. So we rewrite (2.1) as 

S (d) = -\j[(l-5 2 )K(x-y) + 5 2 5(x-y)]F; u (x)U ab (x,y)F^ b (y)dxdy 

= Sk d) + Sl d \ (2.4) 

with 

S ( d) = ~j K(x - y)G; u (x)G^ a (y)dxdy. (2.5) 

Clearly (2.4) is gauge invariant for any 5 and coincides with the original action (2.1) 
for 5 = 1. We used 5 2 and not 5 in the first two terms of (2.4) in order to keep the 
analogy with usual perturbation theory, where the first correction to the wave-function 
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renormalization is also quadratic in the coupling constant. Furthermore the Lagrangian 
(2.4) is symmetric with respect to the transformation — > — A^, 5 — > —8. The free 
action S C ^ is the part of which survives for 5 = 0, while Sj C ^ denotes all the rest. 
Next we perform gauge fixing and quantization with the help of the 
Faddeev-Popov procedure, slightly adopted to our case. We introduce the covariant 
gauge 

F[A^(x)} =fK(x- vWA^y) - X a (y)]dy. (2.6) 
The ghost Lagrangian thus becomes 

sig host) = _ f Q a {x)K{x _ y ){uu a {y) - 6 g f abc d»[A b (y)u c (y)]}dxdy . (2.7) 



At the end of the gauge fixing procedure we integrate over the auxiliary fields x a with 
the weight function 



exp[-^ J X a (x)K(x - y) X a (y)dxdy}. (2.8) 



The result then combines with the terms involving d' 1 d u A® in the usual way. We thus 
arrive at the action S = Sq + Sj with 

So = ~ / A;(x)K(x - y){nA» a (y) + (1/a - l)d^ A a v {y)}dxdy 

uj a (x)K(x - y)Ouj a (y)dxdy, (2.9) 



S T = S\ d) - 5g f abc J d^ a (x)K(x - y)A b ^y)u J c (y)dxdy. (2.10) 
The free gluon propagator will therefore have the form 



DfM=S ab D»M = 5 ab D{q){g^-i q -^), with £ = !-«■ (2.11) 



Here 



DM ^ wrumrv {2A2) 

with K(q) the Fourier transform of K(x). The free ghost propagator is also given by 
D(q). In fact the form factor in the ghost gluon vertex will cancel against the one in 
the ghost propagator in all ghost loops. This is due to the fact that our Faddeev-Popov 
determinant is just the usual one, multiplied with a field-independent term. The form 
(2.6) which leads to (2.7) has, however, the special virtue that the total action given 
below is invariant under the usual BRS transformation. 

We mention that the special case K = l/( with constant ( would lead back to the 
ansatz in which is, however, only useful if one wants to connect renormalized with un- 
renormalized quantities. If one expresses renormalized quantities through renormalized 
parameters, a constant ( drops out. 
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For the explicit calculations we expand the string U (x, y) in powers of 5 and insert 
the series into (2.4). We thus obtain the following expansion of the action, which we will 
need only up to order 5 2 here. 

S = S + 5S^ + 5 2 S {2) + ■ ■ ■ . (2.13) 

So was already given in (2.9), the interaction terms below are classified in an obvious 
notation, where the indices ins, C, S, Gh denote the origin from insertions (the term 
~~ T f[fi( x ~~ y) ~~ K( x — y)]G a llv {x)G^ va {y)dxdy present in (2.4)), commutator, string, or 
ghost interactions respectively. One finds 

S? = -gof abc Jd fl A:(x)K(x-y)A^ b (y)A^(y)dxdy, 

s (i) = 9o fabc J d ^ {x)K{x _ y){x _ y)P A %sx + (1 - s)y)[d»A vc {y) - d» A» c {y)]dxdyds, 
S ( gI = -9of abc j ' d^ a {x)K{x-y)Al{y)^{y)dxdy. (2.14) 

S\l = ~\j d,A:(x)[5(x -y)- K(x - y)WA™{y) - d»A^(y)}dxdy, 

Sj& = - 9 ^f ab J ade J Al(x)Al(x)K(x-y)A^(y)A^(y)dxdy, 

S { cl = 9of abc f ade J A b ^x)Al(x)K(x - y)(x - y) p A d p {sx + (1 - s)y)d»A™(y)dxdyds, 

4? = -|r b 7 a<fe / d,Al(x)K(x - y)Q(s - s')(x - y) p A c p {sx + (1 - s)y) 

(x - y) a A d a {s'x + (1 - s')y)[d^A ve (y) - d v A^{y)\dxdydsds' . (2.15) 

Due to the strings, the action is in fact an infinite series in 5. This is a necessary 
consequence of the non-locality and the gauge invariance of the approach. In any finite 
order we will, of course, only need a finite number of terms. 

The vertices of the action can be visualized by slightly modified Feynman graphs 
which are shown in fig. 1 in momentum space. The propagators refer to D(q) now, the 
presence of a thick line denotes a factor K(q). If such a factor appears in an internal line, 
it results in a total propagator D(q)K(q) = l/(q 2 + ie), i.e. the usual free propagator. 
Thick lines with n gluon lines attached at the interior of the line denote the n-th order 
expansion of the string U(x,y). The graphs make the structure of the action more 
transparent, we did, however, not set up general modified Feynman rules here, but 
preferred the direct one-loop calculation. 
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3 Vacuum polarization 



We start with the one particle irreducible contributions to the full gluon propagator 
which make up the vacuum polarization H^. Due to the nonlocality of So, it is most 
convenient to write down the path integral representation for the propagator, expand it 
with respect to 5 , perform the Gaussian integrations and calculate the contractions in 
momentum space. The whole calculation can now be done using familiar methods. The 
result, written in d dimensions and for N c colors, has the form 

= 5 2 [K(q) - l](gV - <W&0 " "Ap / ^(?, k)d*k. (3.1) 

(2) 

The first contribution at the rhs of (3.1) stems from the insertion S in g in (2.15). The 
integrand n pv (q, k) is a sum of 20 terms which we denote by irffl. Terms (1) to (14) are 

loop terms arising from in (2.14): (l)-(4) from S^S^, (5)-(9) from S^S^, 

(10)-(13) from S^S^, and (14) from S^S^l 

Terms (15)-(20) originate from the second-order part S^ in (2.15): Term (15) is the 
tadpole from the commutator term Sqq, (16), (17) stem from S^s > an d (18) _ (20) from 

(2) 

Sgg . Wherever possible, we simplified the expressions by using relation (2.12). The 
contributions are shown in graphical representation in fig. 2. The explicit forms read: 



k 

i 2 q{q - k) 



k 2 (q — k)'- 



K\q K ~ qkq u ]}D(k)D(q - k)K\q), 



= ^{qkg^ - q^K - £ ^ 9 _ ^ [g(g ~ %^ - q»(q - k) v ]}D(q - k )^r, 

1 



= -(d - l^q - k)„ _ ., 

{q-kf 
K{k) 



= {k 2 9^ + (d- 2)k p k u - — — [(k(q - k)) 2 g^ u - 2k(q - k)q p k v 



+ (q*-k 2 )k ll K}}D(q-ky- 



= {q 2 k p - qkq p }D(k)D(q - k)K{q) J K v {k - sq)ds, 
^jS = Hk(Q- k )9^-q^q-k) u ](q p -^k p )-(q 2 k^-qkq^(g up -^ 1 ^)} 
D(k)D(q - k)K(q) J K p (q - sk)ds, 

7$ = -2{(d - 2)k(q - k)k p + k 2 (q - k) p }^^- J K u (k - sq)ds, 

*J2 = H[q(q-k)g^-qM-kU[(q-k) P -^^^-k p ] 
-(q - k )Mq - k)gv P - (? - k )vq P ] 
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^(14) 



^(16) 



^(18) 



^(19) 



(g - k)' 
D(q-k) 



= 2(9 - fc) M {*?^p - M,} lg fe2 J / ^ P (g - sk)ds, 



Trg?) = i{(rf-2)(%-A;)) 2 + A; 2 (g-A;) 2 } J D(A:) J D(g-A;) 



y y K p (sq — k)K u (s'q — k)dsds\ 

rfu ] = -2{fc(g - fc) [(g - fc)00w> - (g - fcWl + [g% - g 2 ^] (g - k) p } 

D(k)D(q -k) J J K p (q - sk)K u (k - s'q)dsds', 
7r S' 2) = (q - k) u q a }D(k)D(q - k) 



y J K p (q- sk)K a ((l- s')q + s'k)dsds' : 



k,,k n 

,2/ 



+ [g (g - A) M - g(g - k)q^}(q - fc)„)} 
D{k)D(q -k) J y i^(g - sfcjlT^g - s'k)dsds', 

k M - k )v 
k 2 (q — k) 2 ' 



<, 5) = -{(d-i-0^ + e (g ( ^ )2 ^ (g- fc)^(Jb), 



: -2{^g P - - ^(qkg^u - q^k v )k p }D{k) J K p (q + sk)ds, 

- -2(d - l)k^D(k) y K v (k + sq)ds, 

: -[?V ~ qM(9 P « - ^)D(k) J J Q(s - s')K»°{q + (s - s')k)dsds>, 

- -(d - l)k 2 D(k) y y 6(s - sf)Kp,(k + (s - s')q)dsds', 

vrj 20 ) = 2[kqg^-k^]D{k) j j Q{s-s')K pu {(l-s')q + sk)dsds'. (3.2) 

The integrations over s or s' in the string terms, run from to 1. We introduced the short 
hand notations K Pj (k) = d p K(k) and K pu (k) = d' x d v K{k). The string terms contain 
only derivatives of if because they have to vanish for the usual local form if = const. 

Our procedure in the last section guarantees that IL^ is transversal. We check this 
by contracting IL^ with q^q v . We denote the corresponding integrands by tt^, ■ ■ ■ , vr^ 20 ^. 
The non vanishing terms are 

7rf } = -(d- 1)^ -fc) ' 



fc 2 (g — fc) 2 ' 



-; 4) = {g 2 fc 2 +(d-2) ( gfc) 2 - ^ [g ; ~ r j } 



£ 2 [g 2 fc 2 - (gfc) 2 ], L>(g- fc)AT(fc) 



(g - A;) 2 J A; 2 
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D(q — k) 

-2{(d - 2)qk k(q - k) + k 2 q(q - k)}[ 



tt, (7) = -2{(d - 2)qk k{q - k) + k 2 q{q - k)} -^- ® [ q l 'KJk - sq)ds 

k 2 J 

f D( q-k)K(k) 1 

~~ k 2 k 2 {q-k) 2 

ttJ 10) = i{(d-2)(%-fc)) 2 + A; 2 ( g -A;) 2 }^(W?-^) 



y y q^K^sq - k)q v K v (sq - k)dsds' 
{(d-2)(k(q-k)) 2 + k 2 (q-k) 2 }[ 



A; 2 k 2 (q-ky 

(u) _ gfcg(g - k) 
k 2 (q — k) 2 ' 

^ (15) = ~{(d - 1 - O^ 2 + C ^J^ 2 }D( q - k)K(k), 
tt{ 17) = -2(d - l)qkD(k) J q v K v {k + sg)rfs, 



(19) 



k 2 k 2 (q~k) 2i ' 
(d - l)k 2 D(k) J j 6(s - s')q»q v K^(k + (s - s')q)dsds' 

(A 1\/2/ , N2r - P (g- fc)-^(fc) 1 i /q ox 

= -(rf-i)fc (g-fc) [ j^r^l- ( 3 - 3 ) 

The second forms for the string terms 7r} 7 \ 7^°^, 7Tj 17 \ 7Tj 19 ^ which do no longer contain 
integrations over s, s' are obtained as follows: E.g. in n\ 1 ^ use q^q v K^ v (k + (s — s')g) = 
— (d/ds)(d/ds')K(k + (s — s')g), perform the integrations over s and s', and substitute 
k — > /c — p where appropriate. The cancellation of the sum £? 7r F' > now happens in the 
following way: All the ir\^ contain either the factors D(q — k)K(k)/k 2 or the products 
of the free propagators l/k 2 (q — k) 2 . The terms of the first type add up to zero, the sum 
of the latter ones may be written as 

' {d - 2){ Jq-W + ¥ } - (3 ' 4) 

Obviously this expression vanishes after integration over k. 
Due to the transversality of we may write 

IW?) = (Q 2 9^ ~ q^Tliq 2 ). (3.5) 

Il(g 2 ) has the form 

U(q 2 ) = 5 2 [K{q) - 1] - / n(q, k)d d k. (3.6) 

It is easy to obtain the integrands n^(q, k) appearing in H(q 2 ). Since we have already 
checked the transversality for the sum in (3.2), one can simply calculate the trace: n = 



Proceeding further in the usual way one gets for the gluon propagator 

5 ab ^ v {q) = 5 ab {D^(q) + D fip (q)W x (q)D Xu (q) + DUDUD + •••}. (3.7) 
In matrix notation the geometrical series becomes 

(A) = [(D)-' - (II)]- 1 . (3.8) 

This gives 



A,„(q) = A(q)[g,„ - £(q 2 )^}, (3.9) 



with 



A(?) = (," + *)[*(») -n( 9 )] < 310 > 

and a modified g 2 -dependent (,(q 2 ) which is of no interest here. The gluon wave-function 
renormalization constant Z g i uon , defined at the renormalization scale — Q 2 > 0, is ob- 
tained from A(— Q 2 ) = Z g i uon /(—Q 2 ). In order 5 2 this gives 

Z gluon = [1 + Tl{-Q 2 )/K{-Q 2 )]/K{-Q 2 ). (3.11) 

4 Ghost self-energy and vertex function 

A calculation of the renormalized coupling constant from the gluon three-point function 
would be rather complicated. For this one would need the expansion of the action up 
to order S 3 . Furthermore many mixing terms between S^ 1 ) and would show up. 
Therefore we will use the ghost-gluon vertex instead, which is much simpler. 

We start with the ghost self-energy S(g 2 ). It consists only of the graph (21) (fig. 3) 
and reads 

m/f = -^^l* {21 \^k)d% (4.i) 

with 

^>< 9 ,*) = ^-*)-e**V^?S- (4.2) 

It is connected to the ghost propagator A ghost (q) by 

A * w(g) = (<?+ie)[K(q)-Z(q)/q*Y (43) 

Defining the ghost wave-function renormalization constant Z ghost by A ghost (—Q 2 ) = 
Z ghost / (-Q 2 ) one has in order 5 2 

Z ghost = [1 + Z(-Q 2 )/(-Q 2 )K(-Q 2 )]/K(-Q 2 ). (4.4) 
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We next consider the ghost-gluon vertex function T^ c (p, q), where p and q denote the 
momenta of the incoming and outgoing ghost line with color indices a and b. We write 
it in the form 

r?>, q) = / a6 T>, q) = r bc [K(q)q, - J 7|i (p, q, k)d% (4.5) 

There are seven contributions (fig. 4). Graphs (22) - (25) are generalizations of the 
usual vertex graphs, (26) - (28) are string terms. For general momenta the expressions 
are rather lengthy. We therefore specialize to the case p = q, i.e. vanishing gluon 
four momentum. This is sufficient for our purpose. The vertex function simplifies to 
r M (g,g) = q^T(q) with 

T(q) = K(q)- i -^^J 1 (q,k)d d k. (4.6) 

Here 7(9, k) = q^^q, q, k)/q 2 is a sum of seven terms 7^(5, k), only four are different 
from zero. They read: 



(22) r/ M c qkk(q-k)q(q-k)D(k)K(q) 

7 = Ig(g-*)-^fc2— 1 2q^q-kY ' 

^(24) _ r[ 2,2 / k?] q{q-k)D{q-k)K{q) 

7 (23) = 7 (26) = 7 (27) = Q (47) 

The vertex renormalization constant Z vertex is directly related to (4.6): 

Z^rtex = n-Q 2 )- (4-8) 



5 Conceptual questions 

The derivation of the formulae in the last section, though somewhat tedious, was essen- 
tially straightforward. The real problem starts, when the expressions are to be evaluated 
and, in particular, when one has to find a principle for the "optimal" choice of the form 
factor K{k). We will discuss some of the problematics here before entering more detailed 
calculations. 

The evaluation can be performed in the following way. We assume that both K(k) and 
D(k) = l/(k 2 + ie)K(k) satisfy spectral representations with roughly the same behavior 
for large k 2 as in the free case where K(k) = 1, D(k) = l/{k 2 + ie). Therefore we write 
down a once subtracted Kallen-Lehmann representation for K, with the subtraction 
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point chosen at infinity for convenience. For D we assume an unsubtracted dispersion 
relation. Thus put 

K(k>) = 1 + / * {p2) f. , D(k-) = / f^f. . (5.1) 
J k 2 -p 2 + te J k 2 -p? + ie 

Due to (2.12) the spectral functions /t and p are not independent. One could scale K{k) 
with a constant C and -D(A:) with 1/C, respectively. The normalization constant C 
cancels, however, in all the loop graphs contributing to the gluon vacuum polarization, 
the ghost self-energy, and the vertex function. It only enters in the insertion to the 
vacuum polarization (first term on the rhs in (3.1)). In the following the normalization 
constant C will be of no importance, therefore we choose it equal to 1. 

The momentum integrations / • • • d d k in the expressions of the previous sections can 
now be performed in the usual way, although this becomes somewhat ugly for some of 
the string terms. At the end one is left with some integrations over the spectral functions 
R and p. After having calculated, say the vacuum polarization, the ghost self-energy, 
and the ghost-gluon vertex in this way, one has to decide how to choose the "optimal" 
spectral function R, i.e. the input function K(x — y) in our ansatz (2.4). Two well 
known and successful principles suggest themselves: The principle of fastest apparent 
convergence (FAC) postulates that the considered quantity Q does not change when 

going to a higher order of perturbation theory, i.e. Q n = Q n -i- The principle of minimal 
sensitivity (PMS) requires that the quantity be stationary with respect to an arbitrary 

parameter A, i.e. dQ n /d\ = 0. In our case PMS would not simply lead to an extremal 
problem but to a variational problem, because, instead of a single parameter A, we have 
an arbitrary function R at our disposal. 

In a finite theory one could apply FAC or PMS to, say, the propagator and thus obtain 
a non-perturbative solution. In fact this can be done very easily in a toy model like four- 
dimensional $ 3 theory (I thank I. Solovtsov for suggesting this simple exercise and N. 
Brambilla and A. Vairo for an enlightening discussion of the result). FAC leads to an 
integral equation of the Dyson- Schwinger type which can easily be solved by iteration. 
PMS, on the other hand, is not applicable in one-loop order, because the propagator 
becomes independent of the spectral function. 

Fundamental conceptional questions arise, however, if one tries to apply these ideas 
to field theories with divergences as in the present case. It appears natural to apply the 
above principles (FAC or PMS, respectively) to renormalized quantities, expressed by 
renormalized parameters. But the situation is more subtle here than in usual pertur- 
bation theory. The reason is, that renormalized quantities are not necessarily finite for 
a general function K\ This is rather obvious, e.g. from the expressions in (3.2): The 
divergent terms have different g 2 -dependent factors in front and thus will not cancel in 
the differences appearing, say, in renormalized propagators. Technically, the reason is 
that our action (2.9), (2.10), though gauge invariant for every 5, is not renormalizable if 
5 is considered as the coupling constant. 

One way to overcome the problem of divergences would be to simply drop the di- 
vergent contributions in renormalized quantities. From ordinary perturbation theory we 
know that they have to be absent, so also in our case they have to cancel when the whole 
perturbation series in 5 is summed up. The removal of the divergent terms is certainly 
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not unique, analogous to the ambiguities of renormalization schemes. In finite orders 
the results will depend on the detailed prescription, the exact result should, however, be 
independent of it. The remaining finite quantities could then be determined by using 
FAC or PMS. We shall not investigate this possibility here but proceed with a calcula- 
tion of the anomalous dimensions and the /^-function, where the problem of removing 
divergences does not show up. 



6 Anomalous dimensions and /3-function 

The /5-function describes the scaling behavior of the quantized field theory which differs 
from the naive expectations from classical scale invariance. In order to maintain the clas- 
sical scale invariance in every order of the (^-expansion we make a scale-invariant ansatz 
for the spectral functions R (dimension 0) and p (dimension -2). The only scale available 
is the external momentum q in the propagators or the vertex function, respectively. With 
Q 2 = —q 2 > the euclidean squared momentum, we therefore put 

R(jj, 2 ) = K(p 2 /Q 2 ) = K(m 2 ). (6.1) 

The function k as well as the integration variable m 2 = p 2 /Q 2 are now dimensionless. 
From (5.1) we get 

m 2 ) = i + / " (m V • dm2 - ^ 

J k 2 /Q z — m z + xe 

In particular we have 

K(-Q 2 ) = 1 - / ^\dm 2 , (6.3) 
J 1 + m z 

i.e. K(—Q 2 ) becomes independent of Q 2 . Therefore the corresponding factors present 
e.g. in tt^\tt^ 2 \ ■ ■ ■ are not differentiated when calculating the /3-function and the diver- 
gent contributions disappear as usual. 

From simple dimensional analysis in d = 4 — 2e dimensions the integrals appearing 
in (3. 6), (4.1), (4. 6) are proportional to (Q 2 )~ e = 1 — eln Q 2 . Only the divergent terms 
~ 1/e survive after applying the operator Q 2 d/dQ 2 . By naive power counting in k one 
would conclude that all the string contributions in the vacuum polarization (3.2), with 
the exception of -k^ 7 \ tt^ 17 \ should be convergent. This is, however not true. The 
string parameter s has to be integrated from to 1, and for s = most of the terms are 
divergent again by power counting. In fact, a careful investigation allows to extract the 
divergent contribution arising from the region of small s. In the appendix we give some 
technical details, how this can be done in a rather straightforward way. 

The propagators D will be treated by expressing D through K by using (2.12), 
introducing (5.1) (or (6.2)) for K and expanding with respect to the integral: 

dw = w + *) W - (^[i - / + oa/n (6.4) 

The terms of order 1/k 6 lead to finite integrals and are therefore irrelevant. 
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We write the renormalization constants, derived in (3.11), (4.4), (4.8) in the form 



Z 9 luon = ^ 1 + 62 ( K - 1 )/ K+ ^^(- e -^Q 2 )^ 

Z 9host = l[l + *^(I-ln<m (6.6) 

In the above equations we have abbreviated 

K = K(-Q 2 ). (6.8) 

We have only written down the divergent terms ~ (1/e — InQ 2 ) of the integrals 7r, a, 7, 
the factors in front have been denoted by it, a, 7. The renormalized coupling constant 
becomes 



9 = ^Z]ll n Z ghost Z^ rtex = i|[l + ^(K - 1) + 0^(~ e ~ InQ 2 )/?], (6.9) 
with /3 = 7r/2 + o- + 7. In lowest order one has 

g = 5g /VK. (6.10) 
The various contributions to it, a, 7, (5 all have the form 

h + J l(m 2 )K(m 2 )dm 2 + J J K(m 2 ,m f2 )K(m 2 )K(m' 2 )dm 2 dm' 2 , (6.11) 

i.e. they are at most quadratic in the spectral function /t(m 2 ), higher powers only 
contribute finite terms. Below we give the non-vanishing contributions, where the indices 
correspond to the numbers in (3.2), (4.1), and (4.7) as well as to the graphs in the figures. 
All calculations can be easily done analytically, with the exception of the quadratic 
contributions of the string-string graphs (12), (13). We leave them unspecified here. 

It is important to note, however, that graph (13) also has a ^-dependent term which 
can be calculated analytically. This will be important for the cancellation of the in- 
dependence in the /3-function. The non-vanishing contributions read: 









= i+e/2, 




= -1/2, 




= C/2, 




= 1/6, 




= (l + C/2)/2, 


h (22) 


= (1-0/8, 


h (25) 


= 3(1-0/8, 
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(1) 2(1 - e/2) 

1 + m 2 

(2) (l+j/2) 
1 4- W 

Z<-> = - 2m2+ f + 2^1nl±^. 
1 + 

/( 7 ) = 4, 

= - 2m2+ 5 +2^1nl±^. 
1 + m 2 

(g) 1 + 2m 2 2 1 + m 2 

/ (9) = --2m 2 \n — , 

1 + m 2 m 2 

1^ = -4, 

J(18) : 
Z (21) : 
Z (22) : 



" 25 * - "^j- < 6 ' 13 > 



2 + ^ 




1 + m 2 




1 + 




2(14 


- m 2 ) 


1 - 


' 


8(14 


- m 2 ) 


3(1 


-o 



1-^/2 
(l4-m 2 )(l4-m /2 )' 



^(6) ,2m 2 4- £ „ 2] ! + "i 2 1 2 /2 

l + m z m 2 2(1 4- m /2 ) 

#(12) = ^ 12 \m 2 ,m' 2 ), 

K (is) = K^\m 2 ,m' 2 )-— J- -. (6.14) 

v ' ; 2(1 4-m 2 )(l 4-m' 2 ) v ; 

Summing up contributions (1) - (20) for rr, (21) for a, and (22) - (25) for 7, we obtain 

h , = 5/3 + { /2, i , = -5^2 + 2^h 1 i±^ 

1 + m 2 m 2 

1 4- m 2 4- m' 2 m 2 1 + m 2 m' 2 1 4- m' 2 
K n = — ^rrz — — ~ In m 



(1 4- m 2 )(l 4- m' 2 ) 1 + m' 2 m? 1 + m 2 m' 2 
+K( 12 \m 2 ,m> 2 ) + K^ 3 \m 2 ,m' 2 ), 



(6.15) 



^ = (l + e/2)/2, l* = - 2 \i + £ 2 y ^ = > ( 6 - 16 ) 
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^ = (1-0/2, ^ = - 2(1 1 + ^ 2) » K ^ = °- ( 6 - 17 ) 



Furthermore we have for (3 = n /2 + a + 7: 



1 + m 2 

The anomalous dimensions and the /3-function become 



hp = 11/6, lp = m 2 \n—^-l, K p = -K^. (6.18) 



' r\2 gluon 
7 gluon — V Apil 

^gluon OiUf 



g 2 N c r 5 C, f , , 2 \ , 2\ 1 2 
+ - + y /^(m )/t(m Jam 



(4tt) 2L 3 

+ / [ K n (m 2 ,m' 2 )K{m 2 )K(m' 2 )dm 2 dm' 2 }, (6.19) 



*~ - Zb ^ - "ffi + ! + /«™>(™ 2 )M, (6.20) 



-'ghost 



dQ 2 (4tt) 



^ 2 « 2 ^ ^ -^[T + /W» a )«< m ')*» a 

+ | [ Kp(m 2 ,m' 2 )K(m 2 )K(m' 2 )dm 2 dm' 2 }. (6.21) 

We used (6.10) to replace <5go/ by the renormalized coupling constant g. For /t(m 2 ) = 
we recover the well known results of ordinary perturbation theory in one-loop order. 
Note further that the /3-function is independent of the gauge parameter £ as it should. 
This is a further test of the manifest gauge invariance of the approach. 

Let us next discuss whether the principle of fastest apparent convergence (FAC) or the 
principle of minimal sensitivity (PMS) can be applied, say, to (3 in (6.21). An application 
of FAC is clearly impossible because the present one-loop calculation is the lowest non- 
trivial contribution; there is no lower order with which one could compare. An attempt 
to apply PMS also fails. The variation with respect to K,(m 2 ) gives no solution at all in 
(6.20) which is linear in k. In (6. 19), (6. 21), on the other hand, the presence of the linear 
term / l(m 2 ) K{m 2 )dm 2 leads to an extremum which is not situated at k — 0. Therefore 
one would not reproduce the results of ordinary perturbation theory for small g. 

This is not yet a problem. We know from simple toy models that the lowest order 
usually gives no extremum. One needs the lowest and, at least, the next order to get a 
balance between these contributions and to find a relevant extremum. In both cases, FAC 
or PMS, one should therefore go to order S 4 for the anomalous dimensions and to order S 5 
for the /3-function. A calculation of the two loop contributions ~ 5 5 g 5 to the /3-function 
is not feasible at present, but one can easily look for the effect of the contributions ~ 5 5 g 3 
which, of course, are the most important ones for small coupling. These contributions 
arise from two sources: First we have to consider the insertion 5 2 (K — 1)/K in Z g i uon , 
when expressing 5g by g. This leads to 
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9 = ^=[1 + y (1 " V*) + (6.22) 



Inverting this equation, one obtains for an arbitrary power n 



%i =' + + (6.23) 



\VKj a 2 
Thus, for 5 — 1, instead of (6.10), one should now replace 

Different replacement rules for different powers look a bit strange but appear as a direct 
consequence of the expansion in 5. 

The second effect is that we now have to consider insertions into the internal gluon 
lines in the graphs (l)-(28). Such an insertion replaces the original propagator D^ u {k) = 
D(k 2 )(g^-S > k lJL k u /k 2 ) by (1 - 1/ K(k 2 ))D(k 2 )(g lJLU - k^k u /k 2 ). This is transversal, there- 
fore it is most easily discussed in the Landau gauge £ = 1 where one simply gets a factor 
1 — 1/K{k 2 ) = J dfi 2 R,(fi 2 )/(k 2 — fj, 2 + ie) + 0(l/k 4: ) which multiplies the original propaga- 
tor. Obviously one has a suppression by an additional power of two in the denominator, 
therefore only insertions into the graphs (3), (4), (15) lead to divergent contributions, all 
the others become convergent. In Landau gauge we find an additional contribution of -4 
to l (3 \ -4 to 1 {A \ and 12 to /( 15 ). 

An immediate result is, that now ^ghost becomes independent of the spectral function 
K(m 2 )\ because of the relation l a = —h a /(l + m 2 ) in (6.16) the square bracket in (6.20) 
is proportional to 1 — / dm 2 K(m 2 ) / (1 + m 2 ) = K which cancels against the 1/K in 
front which now survives when 5 2 gl/K is replaced by g 2 / K according to (6.24). This 
is quite welcome, because the unpleasant linear contribution in k has thus disappeared. 
The same is, however, not true for 7 9 ; uori and for (5. This result will also hold in higher 
orders, because higher order insertions make all graphs finite. 

The additional contributions of order S 5 g 3 just discussed (remember that we did not 
take into account terms of order 5 5 g 5 ) all vanish if K,(m 2 ) = 0. Comparing order 5 3 
and <5 5 , FAC would now trivially give the solution k = which is just what we expect. 
Inclusion of the two-loop contributions ~ 5 5 g 5 would shift this to a non-trivial solution 
for k and finally result in a non-perturbative solution for the /3-function. 

For PMS, on the other hand, there is really a problem. The contributions ~ k are still 
there, and the optimized (^-expansion would not reproduce the lowest order perturbative 
result! 



7 Conclusions 

At present we cannot offer a convincing explanation for the results obtained above. Of 
course we cannot exclude the possibility of a calculational error in our formulae, although 
we have carefully checked them. Assuming that they are correct, it remains a puzzle 
why PMS fails to reproduce the lowest order perturbative result. This is particularly 
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confusing because the puzzle persists to any order in 5 if g is small. An encouraging 
result would have consisted in a cancellation of all the terms linear in K(m 2 ) appearing 
in the /3-function, i.e. lp = 0. The quadratic terms would then lead to an extremum at 
K(m 2 ) = 0, thus reproducing the results of one-loop ordinary perturbation theory. In 
a two-loop calculation, which might be feasible with more intensive computer help, one 
would then expect a non-trivial solution for K,(m 2 ) and a non-perturbative determination 
of the /5-function. But this was not what we obtained. On the other hand, FAC could 
lead to interesting non-perturbative results in two-loop order which include ordinary 
perturbation theory when expanded with respect to g. 

Let us now discuss a different possible approach which was already briefly mentioned 
in sect. 5. One could calculate the renormalized gluon propagator in terms of the general 
spectral function R introduced in (5.1) and remove the remaining divergent contributions 
by a definite prescription. This would be complementary to - and of course much more 
complicated than - the calculations in sect. 6 where we concentrated on the divergent 
terms. Finally one could determine R using FAC or PMS. 

Note, however, that our ansatz in (5.1) implies K(k 2 ) — > 1 and, accordingly, D(k 2 ) — > 
1/k 2 for k 2 — > oo. The renormalized propagator, as determined by PMS, could neverthe- 
less show the correct asymptotic behavior as obtained from the renormalization group, 
namely A(k 2 ) ~ [In k 2 }-^+^' AA /k 2 . (See e.g. §). If one would apply FAC, on the 
other hand, the bare and the renormalized propagator would be identical by definition 
of the method. In this case one should either use the special gauge £ = —10/3 for which 
Igiuon = and A(A; 2 ) ~ 1/k 2 , or, more generally, write down a twice subtracted disper- 
sion relation for K(k 2 ) with a finite subtraction point. A calculation of this type, though 
complicated, appears possible and will be undertaken in the future. 

Finally we would like to mention a conceptional problem. The connection between 
the bare and the renormalized coupling constant will always start with g ~ 5 go . When g 
is eliminated, higher powers {Sgo) n become proportional to g n , i.e. become independent 
of S. So the clear bookkeeping of powers of S is somehow blurred by the renormalization 
procedure. 

The present approach appears complicated. We interpret this as a reflection of the 
fact that QCD is complicated and that non-perturbative results can only be obtained 
with considerable effort. We believe that interesting information can be extracted from 
the general expressions presented here, and even more interesting information from a 
two-loop calculation which might be feasible. The puzzle in connection with PMS is not 
understood at present. Any suggestions are welcome. 

Acknowledgement: I thank N. Brambilla, M. Jamin, G. Prosperi, I. Solovtsov, 
and A. Vairo for valuable discussions and I. Bender for a careful reading of the manuscript. 
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A Appendix 



We give here some technicalities how the divergent contributions of the string terms 
can be extracted. The problematic terms are those where the internal momentum k is 
multiplied with a string parameter s which has to be integrated from to 1. For s ~ 
there is no suppression by a power of k in the denominator, therefore a more detailed 
analysis is necessary. 

In a first step introduce the spectral representations (5.1) or (6.2), (6.4) for K and D, 
as well as the resulting representations for and K^, and rotate to euclidean space 
as usual. In all denominators without string parameters s or s' expand with respect to 
1/k 2 , i.e. 

1 1_ 2qk 4(qk) 2 q 2 + /i 2 

(k - q) 2 + /i 2 k 2 k 4 k 6 k 4 1 ' ' 

up to the order where the further terms become ultraviolet finite. The apparent infrared 
singularities arising from the expansion are spurious. 

In terms with only one string, e.g. (8), we have a further denominator of the form 
l/[(q—sk) 2 +fi 2 ] 2 arising from the spectral representation of K p (q — sk). The substitution 
k = k'/s makes the integration over k' finite. But for every power k~ n in the integrand a 
power s n ~ d = s n ~ 4+2e appears in d = 4 — 2e dimensions. The s-integration gives a factor 
(n — 3 + 2e) _1 and thus a divergent contribution if n = 3. 

Next we discuss the terms with two string integrations over s, s'. Terms (10), (11) are 
finite, in (12), (13) substitute s = r(l—t), s' = rt. This gives a factor r from the Jacobian. 
Next substitute k = k'/r and perform the integral over r from to 1 as before. This gives 
again some divergent factors, the remaining integral is finite. The second order term in 
(19) is finite, in (18) and (20) substitute t = s - s' ,t' = (s + s' - l)/2. The t'-integration 
can be trivially performed and gives a factor (1 — t), then substitute k = k'/t as before. 
After integration over k! the various terms in (20) cancel, so that this graph, contrary 
to what one would expect, does not give a divergent contribution. 
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Figure Captions 

Fig. 1: Interaction terms in momentum space. All momenta are incoming. A thick 
line denotes the presence of a factor K. Thick lines with gluons attached at the interior 
of the lines arise from the expansion of the string U (x, y) . They are associated with a 
factor = d^K or K^ v = d tx d v K, respectively. String parameters s, s' are integrated 
from to 1. 

Fig. 2: Contributions to the gluon vacuum polarization. 
Fig. 3: Ghost self-energy. 

Fig. 4: Contributions to the ghost gluon vertex. 
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Fig.2a 
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